Chaos around the superposition of a monopole and a thick disk 
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Abstract 

We extend recent investigations on the integrability of oblique orbits of test 
particles under the gravitational field corresponding to the superposition of 
an infinitesimally thin disk and a monopole to the more realistic case, for 
astrophysical purposes, of a thick disk. Exhaustive numerical analyses were 
performed and the robustness of the recent results is confirmed. We also found 
that, for smooth distributions of matter, the disk thickness can attenuate the 
chaotic behavior of the bounded oblique orbits. Perturbations leading to the 
breakdown of the reflection symmetry about the equatorial plane, neverthe- 
less, may enhance significantly the chaotic behavior, in agreement with recent 

studies on oblate models. 
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The recent observational evidences suggesting that huge black-holes might inhabit the 
center of many active galaxies have motivated some investigations on the dynamics of 
test particles in gravitational systems consisting in the superposition of monopoles and disks. 
Infinitesimally thin disks are frequently used to model flattened galaxies [0. Some exact 
relativistic solutions describing the superposition of non-rotating black-holes and static thin 
disk, and their respective Newtonian limits, were presented and discussed in In [Q, the 
integrability of oblique orbits of test particles around the exact superposition of a black-hole 
and an infinitesimally static thin disk was considered. Bounded zones of chaotic behavior 
were found for both the relativistic and Newtonian limits. There are several examples in the 
literature of chaotic motion involving black-holes: in the fixed two centers problem 0-0] , in 
a black-hole surrounded by gravitational waves MM, and in several core-shell models with 



relevance to the description of galaxies (see [T^ for a recent review). As to the Newtonian 
case, we notice, for instance, the recent work of C. Chicone, B. Mashhoon, and D. G. Retzloff 



on the chaotic behavior of the Hill system |11 



The Newtonian analysis of [^] has revealed an interesting property of the dynamics 
of oblique bounded orbits around the superposition of a monopole and an infinitesimally 
thin disk. Since one was manly interested in bounded motions close to the monopole, it 
was assumed that the disk was infinite and homogeneous. This situation corresponds the 
simplest superposition of a monopole and a disk. Using cylindrical coordinates (r, 6, z) with 
the monopole, with mass M, located in the origin and the disk corresponding to the plane 
2; = 0, the gravitational potential is given by 

M 

V{r,e,z) = -^=== + a\zl (1) 

where a is a positive parameter standing for the superficial mass density of the disk. The 
angular momentum L in the z direction is conserved, and we can easily reduce the three- 
dimensional original problem to a two-dimensional one in the coordinates (r, z) with the 
Hamiltonian given by 



The Hamiltonian is smooth everywhere except on the plane z = 0. Moreover, the parts 
of the trajectories restricted to the regions z > and z < are integrable The corre- 
sponding Hamihon-Jacobi equations restricted to these regions can be properly separated 
in parabolic coordinates, leading, respectively, to the second constants of motion 

Cz>o = Rz-a— (3) 

and 

a<o = Rz + a—, (4) 

where Rz is the z component of the Laplace- Runge-Lenz vector. Note that C is not smoothly 
defined on the disk. With the two constant of motions H and C, the equations for the tra- 
jectories of test particles, restricted to the regions z > and z < 0, can be properly reduced 



to quadratures in parabolic coordinates [0,0- ^ complete bounded obhque trajectory. 



therefore, corresponds to the matching of an infinite number of integrable trajectory pieces. 
Hence, the widespread zones of chaotic motion detected in have their origin in the changes 
in the value of the constant C when the test particle crosses the disk z = 0. As to the rel- 
ativistic case, in contrast, the trajectory pieces that do not cross the disk are themselves 
non-integrable, leading to typically larger chaotic regions than in the corresponding Newto- 
nian limit 1^. 

Here, we study the robustness of the results for the Newtonian case by considering the 
more realistic case, for the description of flattened galaxies, of a superposition of a central 
monopole and a smooth thick disk with the potential (|l|) as the vanishing disk thickness 
limit. A smooth distribution of matter is considered as a disk if its radial gradients are much 
smaller than its vertical ones. A minimally realistic model for a rotating thick disk should 
obey Emden's equation for the stability of rotating polytropes 0. As in 0], we will neglect 
the radial gradients, and, in this case, Emden's equation for the disk matter density p{z) 
states that {G = 1) 

nXp^'-^p" + kX{X - 2)p^'^{p'f = -Awp, (5) 



where k is the parameter relating the pressure to the matter density in the polytropic 
equation of state and A = l + l/ra, n being the polytrope index. The matter density p is 
assumed to obey Poisson's equation V^Vd = 47rp. For the isothermal case (A = 1), equation 
(H) admits as a solution the following distribution of matter 



a , 9 z 



p{z) = sech — , (6) 

Attzq Zo 



which corresponds to the potential 



z 



Vt){z) = 0:2:0 In cosh — , (7) 

where zq measures the disk "thickness" , and a its "superficial" mass density. They obey the 
relation 2k = az^. Typically, for realistic models of a rotating dust disk one has Zq oc {V^) 



14^ . The matter density (y) corresponds, therefore, to a stable and smooth distribution 
of rotating matter concentrated on the plane z = 0. Moreover, in the limit of zq ^ 0, we 
recover from (|^) the infinitesimal thin disk potential Vj){z) = a\z\ and the corresponding 6 
distribution of matter from (^. 

Thus, the dynamics of test particle moving around the superposition of our smooth thick 
disk and a monopole will be governed by the following smooth Hamiltonian 

^2 ^2 jJ2 jy^ ^ 

H = — + — + — ~ ^=== + azolncosh—. (8) 
2 2 2r^ + z^ Zq 

We wish to stress that the potential in (§) corresponds, indeed, to a first approximation of a 
realistic superposition of a monopole and a smooth thick disk with matter distribution given 
by d^). The whole superposition must also obey Emden's equation; in the present case it 
fails to obey at the origin. We are neglecting the radial gradients of the matter distributions 
caused by the stresses induced by the central monopole. 

The Hamiltonian (M) does not belong to the class of integrable two-dimensional potentials 



with a second constant of motion polynomial in the momenta ||T5[. We will present strong 
evidences that (|^) does not have a second constant of motion at all. We could solve numer- 
ically the system governed by (H) with great accuracy. Figure 1 shows a typical Poincare's 



section {H = — 0.2, L = M = l,a = 0.1, 2:0 = 1-5) across the plane z = revealing a 
widespread chaotic behavior. The disk thickness zq, in this case, has the same magnitude of 
the typical z-amplitude of the trajectories. Figure 2 shows a sequence of low-energy sections 
{H = -0.3, L = M = l,a = 0.1, = 0.0 (a), 0.1 (b),0.25 (c), and 0.5 (d)), constructed 
from the same trajectory initial conditions, where the attenuation of the chaotic behavior 
due to the disk thickness can be clearly appreciated. We could obtain thousands of intersec- 
tions for each trajectory with a cumulative error, measured by the constant H, inferior to 
10^^^. We notice that the equation of motion are invariant under the following rescalings: 

r ^ Ar, z^ Xz, A^/^t, 

a — > A~^a, M M, zq Xzq 

H X-^H, L X^/^L- (9) 

and 

r AV, z X'z, t X't, 

a X'~^a, M A'M, zq X'zq 

H ^H, L-^ X'L; (10) 

A > and A' > 0, leading that, for each triple of nonzero parameters (M, a, zo), one has, in 



fact, only one free parameter, namely Zoya/M. 

The limit of zq much larger than the typical z- amplitude of the trajectories deserves a 
special attention. For this case, the Hamiltonian can be well approximated by 

M /3 2 

^- = y + T + 2;:^-7F^+2^' ^''^ 

where P = a/zQ. Such a potential is related to the potential of the Kepler problem perturbed 
by a quadrupole halo potential considered in |T^. Figure 3 shows a typical Poincare's section 



{Hoc = —0.15, L = M = 1, P = 0.1) across the plane z = revealing a widespread chaotic 
behavior for the system governed by (|TT]). We could also obtain thousands of intersections 
for each trajectory with a cumulative error inferior to 10~^^. Like in the infinitesimally thin 



disk case [Q, due to the existence of two rescaling invariances, this Poincare's section can 
be obtained for any non-zero values of (3 and M. 

Our numerical finds confirm the robustness of the results presented in The chaotic 
behavior of bounded orbits can be considered as inherent to any system consisting in the 
superposition of a disk and a central monopole. We stress that the superpositions we have 
considered are symmetric under the reflection about the equatorial plane. A perturbation 
leading to the breakdown of the reflection symmetry could increase significantly the chaotic 



behavior, as the recent studies on oblate models have suggested ||T^ . We could indeed check 



such fact by considering a dipole-like perturbation of the Hamiltonian (H) 

H = H-Dj-^^^. (12) 

Figure 4 presents a typical Poincare's surface {H = —0.3,L = M = l,a = 0.1, 2:0 = 
0.25, D = 5 X 10^^) for this case. The zones of chaotic motion are larger than in the 
corresponding case for which D = (Fig 2.c). However, the central family of regular orbits, 
if it does exit for D = 0, seems to be robust against dipole-like perturbations. Analogous 
conclusions hold also for the large Zq case. 
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FIGURES 



FIG. 1. Poincare's section (r, r) across the plane z = for oblique orbits, with H = —0.2 and 
L = 1, around the superposition of a monopole with mass M = 1 and a smooth disk with surface 
mass density a = 0.1 and thickness zq = 1.5. 
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FIG. 2. Sequence of Poincare's section (r, r) across the plane 2; = for oblique orbits, with 
H = —0.3 and L = 1, around the superposition of a monopole with mass M = 1 and a smooth 
disk with mass surface density a = 0.1 and thickness zq = 0.0 (a), zq = 0.1 (b), zq = 0.25 (c), and 
zo = 0.5 (d). The sections were constructed from the same trajectory initial conditions. 

FIG. 3. Poincare's section (r, r) across the plane 2: = for oblique orbits, with H^o = —0.15 
and L = 1, around the superposition of a monopole with mass M = 1 and a smooth disk of large 
thickness with /? = 0.1. 
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FIG. 4. Poincare's section (r, f) across the plane z = for oblique orbits, with H = —0.3 and 
L = 1, around the superposition of a central source characterized by a monopole charge M = 1 
and a dipole strength D = 5 x 10~^, and a smooth disk with mass surface density a = 0.1 and 
thickness zo = 0.25. The case corresponding to = is presented in Fig. 2.c. 
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